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I. 


INTRODUCTORY  REMARKS 


A  difficult  and  vexing  problem  in  electromagnetic  theory  is  the 
discrimination  of  obstacle  shapes  by  far-field  scattering  measurements. 

The  voluminous  data  usually  taken,  the  ordering  and  reduction  of  these 
data,  the  compromises  which  must  be  made  in  choice  of  measurement 
parameters,  such  as  frequency,  polarization,  etc.,  all  compound  the 
problem  and  make  it  quite  difficult  to  extract  the  shape-dependent 
essence  upon  which  to  base  a  discrimination  scheme. 

Although  this  report  hardly  solves  the  problem  it  does  present 
a  view  point  which  may  be  useful.  This  viewpoint  recognizes  that  a 
good  deal  of  information  about  an  obstacle  is  contained  in  low  fre- 
quencies[l ,2 ,3] ,  i.e.,  those  where  the  obstacle  is  only  a  fraction  of 
a  wavelength  in  maximum  dimension.  At  such  low  frequencies,  it  is 
assumed  that  an  exciting  plane  wave  induces  in  the  obstacle  only  six 
modes  causing  scatter  -  three  electric  dipoles  criented  along  the 
three  orthogonal  axes  of  a  Cartesian  coordinate  system  chosen  within 
the  body,  and  three  magnetic  dipoles  oriented  along  those  same  axes. 

The  relative  strengths  of  these  six  induced  dipoles  are  proportional 
to  the  incident  electric  or  magnetic  fields  colinear  with  them  and  the 
proportionality  constants  may  be  arranged  to  form  a  polarizability  matrix. 
The  terminology,  "polarizability  matrix",  is  familiar  in  the  theory  of 
electrostatic  and  magnetic  static  fields[4,5]  and  has  been  applied  to 
quasistatic  fields  as  well [6,7].  We  have  carried  it  over  to  the 
regime  of  quasistatic  fields  by  relaxing  two  constraints.  First,  and 
in  contrast  with  the  static  case  where  the  polarizability  matrix  has  all 
real  elements,  we  allow  the  elements  of  our  matrix  to  become  complex, 
thereby  permitting  depolarization  effects  other  than  linear-to-linear. 
Second,  and  again  in  contrast  with  the  static  case  where  electric  and 
magnetic  dipoles  are  all  uncoupled,  we  allow  such  coupling  to  manifest 
itself  in  non-zero  off-diagonal  matrix  elements. 

The  36  elements  of  the  6x6  polarizability  matrix  are  determined 
by  illuminating  the  obstacle  in  turn  with  six  different  excitations  - 
each  one  designed  to  most  strongly  induce  only  one  of  the  dipoles 
mentioned  above.  The  far-scattered  field  resulting  from  each  such 
excitation  is  then  calculated  by  computerized  boundary  solution 
technique,  such  as  reaction  matching  introduced  by  Richmond[8] .  The 
reactions  of  each  of  the  six  dipole  fields  with  each  of  the  six 
scattered  fields  (36  reactions  in  all)  are  then  evaluated  over  the  sphere 
at  infinity,  thereby  determining  the  36  elements  of  the  polarizability 
matrix.  (Actually,  matrix  symmetry  reduces  the  maximum  number  of 
necessary  reactions  to  21.) 

Assuming  the  validity  of  the  model,  knowledge  of  the  polarizability 
matrix  permits  the  evaluation  of  the  far-field  scattered  in  any  direction 
due  to  a  plane  wave  of  any  polarization  incident  from  any  direction,  all 
at  a  chosen  frequency.  Such  generality  gives  us  some  confidence  that  the 
polarizability  matrix  of  a  given  obstacle  is  an  almost  unique 


mathematical  representation  of  that  obstacle  ("almost",  because  absolute 
uniqueness  requires  knowledge  of  the  polarizability  matrix  over  a  range 
of  frequencies  -  this  has  not  been  investigated).  If  such  a  one-to-one 
correspondence  between  an  obstacle  shape  and  a  matrix  can  be  assumed,  we 
have  here  a  mathematical  model  which  is  suited  to  processing  by  com¬ 
puter  for  shape  discrimination  purposes.  For  example,  if  [oq]  and  [o^] 
are  polarizability  matrices  representing  two  obstacles,  a  decision¬ 
making  process  could  be  based  on  the  complex  parameter. 


> 


where  [  ]ct  means  "complex  conjugate  transpose  of"  and  '[  ;  means 
"norm  of".  Or,  perhaps  the  elements  of  [a]]  and  [c^]  can  locate  two 
points  in  a  complex  21 -dimensional  space,  the  distance  between  them, 
defined  by  a  suitably  chosen  norm,  serving  as  the  basis  of  a  decision¬ 
making  process. 
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The  last  suggestion  brings  to  fore  a  basic  weakness  of  the 
polarizability  matrix  description  of  an  obstacle,  and  that  is  its 
dependence  on  the  location  and  orientation  of  the  coordinate  system 
chosen  within  it.  Thus,  an  obstacle  is  represented  not  by  a  unique 
matrix,  but  by  a  set  of  matrices  (ideally)  related  to  each  other  by 
rotation  and  translation  operations.  The  best  one  can  hope  for  is 
to  choose  coordinates  which  cause  the  six  dipoles  to  be  characteristic 
modes[9,10,l 1 ,12]  of  the  obstacle,  that  is,  which  cause  the  six  dipoles 
to  be  located  and  oriented  in  such  a  manner  that  each  one  is  excited  to 
a  strength  proportional  to  its  own  normalized  radiation  pattern.  If 
this  can  be  achieved,  then  only  six  complex  numbers  (which  are  unique  to 
within  the  accuracy  of  the  dipole  approximation  in  the  form  of  a 
diagonal  matrix)  suffice  to  describe  the  obstacle  and  can  therefore  be  used 
as  a  basis  for  discrimination.  The  relationship  between  the  polarizability 
matrix  and  the  characteristic  matrix  representations  has  not  been 
developed  sg  far. 

In  this  report  we  present  a  set  of  sample  data  in  graphical  form 
comparing  the  backscatteririg  cross  sections  of  a  variety  of  perfectly 
conducting  wire  obstacles  as  obtained  by  reaction  matching  and  approxi¬ 
mated  by  the  polarizability  matrix  formulation.  The  theoretical 
development  of  this  formulation  is  given  in  the  Appendix. 
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II. 


NUMERICAL  RESULTS 


Numerical  results,  comparing  backscattering  patterns  calculated 
using  Richmond's  reaction  matching  technique  and  the  polarizability 
matrix,  were  calculated  for  several  simple  configurations  -  the  linear 
wire,  the  linear  wire  tilted,  the  square  loop,  a  cube  and  a  rec¬ 
tangular  box,  both  outlined  by  wires. 

figures  1-3  show  backscattering  cross  sections  of  three  lengths  of 


Fig.  1.  Backscattering  cross  section  as  a  function  of 
incident  angle  e  for  a  dipole  of  length  =  0.1A 
and  radius  =  0.005a.  The  incident  electric 
field  is  9-polarized. 
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Backscattering  cross  section  as  a  function  of 
incident  angle  a  for  a  dipole  of  length  =  0.2 
and  radius  =  0.005A.  The  incident  electric 
field  is  6-polarized. 


Fig.  3.  Backscattering  cross  section  as  a  function  of 
incident  angle  e  for  a  dipole  of  length  =  0.5; 
and  radius  =  0.055A.  The  incident  electric 
field  is  8-polarized. 
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9 ( OEGREES ) 

Fig.  5.  Backscattering  cross-section  as  a  function  of 
incident  angle  0  for  a  tilted  dipole  of 
length  =  0.1A  and  radius  =  G.005A.  The  inci¬ 
dent  electric  field  is  <t>-polarized. 

data  for  a  tilted  wire  of  length  A/10  illuminated  in  two  polarizations. 
Similar  data  for  a  wire  of  length  A/5  are  shown  in  Figs.  6-7. 

Figures  8-9  show  backscattering  patterns  of  a  square  loop,  A/10 
on  a  side  in  two  polarizations,  and  Figs.  10  and  11  show  similar  da*-3 
for  a  cube  A/25  on  a  side  and  a  rectangular  box  A/25  x  A/25  x  A/15  in 
dimension.  Both  of  these  latter  obstacles  are  not  solid  bodies  but  wire 
outlines  as  indicated  in  the  figures. 

All  cases  presented  show  good  agreement  with  the  "exact"  solution 
of  Richrr  id.  Although  not  plotted,  phase  patterns  are  equally  good. 
However,  for  obstacles  larger  than  the  ones  shown,  the  polarizability 
matrix  representation  quickly  loses  accuracy,  particularly  for  the 
more  complex  shapes.  The  limitation  to  cubes,  f 0  *  example,  about  A/25 
on  a  side  is  disappointingly  restrictive  and  severely  limits  the 
practical  use  of  the  representation. 
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Fig.  6.  Backscattering  cross  section  as  a  function  of 
incident  angle  6  for  a  tilted  dipole  of 
length  =  0.2A  and  radius  =  0.005X.  The  inci¬ 
dent  electric  field  is  6-polarized. 
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Backscattering  cross  section  as  a  function  of 
incident  angle  8  for  a  tilted  dipole  of 
lengtn  =  0.2X  and  radius  =  0.005X.  The  inci¬ 
dent  electric  field  is  ^-polarized. 
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Fig.  8.  Backscattering  cross  section  as  a  function  of 
incident  angle  8  for  a  square  loop  of 
edge  =  0.1  X  and  wire  radius  =  0.0Q5A.  The 
incident  electric  field  is  8-polarized. 
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Backscattering  cross  section  as  a  function  of 
incident  angle  6  for  a  square  loop  of 
edge  =  0.1  A  and  wire  radius  =  0.005 A.  The 
incident  electric  field  is  ^-polarized. 
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Fig.  10.  Backscattering  cross  section  as  a  function  of 
incident  angle  e  for  a  cube  cf  edae  =  0.04X. 

(a)  The  incident  electric  field  is  e-polarized. 

(b)  The  incident  electric  fie^d  is  $ -polarized. 


III.  CONCLUSIONS 


For  electrically  small  obstacles  (by  which  we  mean  obstacles  at 
most  a  fraction  of  a  wavelength,  i.e.,  A/5  -  A/25  in  maximum  dimension) 
the  polarizability  matrix  representation  is  an  accurate  description  of 
the  scattering  properties  of  obstacles.  This  size  limitation  is  smaller 
than  envisioned  for  obstacles  pertinent  to  the  objectives  of  the  present 
grant.  The  addition  of  multi  Doles  of  higher  order  than  the  dipole 
would  improve  the  representation  and  probably  would  extend  its  accuracy 
to  larger  boa5es,  but  at  severe  expense  of  simplicity.  Additionally, 
there  remains  the  question  of  how  discriminatory  is  the  polarizability 
matrix  representation,  that  is,  how  "different"  are  the  matrices 
associated  with  obstacles  of  differing  shapes,  so  that  they  can  be 
used  as  sensitive  and  reliable  discriminators  of  shape.  Further 
investigation  into  such  questions  is  of  interest,  but  cannot  be 
recommended  at  the  present  time  in  view  of  the  aims  of  the  present  grant. 
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APPENDIX  I 

THEORETICAL  DEVELOPMENT 


At  very  low  frequencies,  the  fields  scattered  by  an  obstacle 
arise  from  rwc  induced  dipole  moments,  one  electric  and  one  magnetic, 
each  of  which  may  be  decomposed  into  three  dipoles  along  three  mutually 
perpendicular  axes  of  a  Cartesian  coordinate  system  (Fig.  1-1)  choser. 


z 


Fig.  1-1.  Polar  coordinate  system. 


within  the  obstacle.  The  six  induced  dipole  moments,  pex,  pey,  pe2, 
Pmx»  Pmyj  Pmz>  (^9-  I-2)  are  complex  and  are  related  to  the  exciting 
electric  and  magnetic  fields  (E_i  .H1 )  by 


where  (F1 )  = 


X 

E1 
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ri 
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is  a  column  vector  composed  of  the  Cartesian 
components  of  the  exciting  electric  and  mag¬ 
netic  fields  at  the  chosen  origin.  They  may 
be  complex. 


[a]  is  defined  to  be  the  polarizability  matrix. 


If  the  electric  and  magnetic  fields  are  those  of  a  plane  wave,  then 
the  vector  (P)  is  understood  to  depend  upon  the  polarization  state 
and  direction  of  incidence  of  the  plane  wave  with  respect  to  the 
chosen  coordinate  system.  The  components  of  the  polarizability  matrix 
depend  only  upon  the  obstacle  shape  and  composition  and  the  coordinate 
choice. 

In  matrix  notation,  the  far  electric  field  pattern  evaluated  in 
the  (es,$s)  direction  may  be  obtained  from 
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[t  =  [R]  (p). 
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where  Z0  =  120tt  and  the  constant.  g~  normali zes  the  R  matrix  such  that 

a  unit  induced  dipole  moment,  e.g„,  pex,  radiates  unit  power.  Combining 
Eqs.  (2)  and  (4),  we  obtain  an  expression  for  the  far  field  scattering 
pattern  due  to  a  plane  wave  arriving  from  any  direction, 

(5)  (E5)  -  [R]  f«]  (f1 ). 


To  determine  the  elements  of  the  polarizability  matrix,  we 
illuminate  the  obstacle  with  6  different  excitations,  each  of  which 
causes  the  (F1)  matrix  to  exhibit  only  one  non-zero  component.  Each 
such  excitation  consists  of  two  plane  waves  traveling  in  opposite 
directions  along  the  3  principal  coordinate  directions,  and  polarized 
constructively  or  destructively  to  produce  a  null  in  either  the  electric 
field  or  the  magnetic  field  at  the  origin.  The  six  excitations  are 
pictured  in  Fig.  1-3.  The  far  electric  field  patterns  scattered  by 
each  excitation  can  be  calculated  by  an  independent  method  such  as 
Richmond's  reaction  matching[8]  and  are  symbolized  by  EgX,  E^y,  EgZ, 

EmX,  Ejny,  Ejjjz,  where  the  subscripts  denote  tne  type  and  direction  of 
the  field  produced  at  the  origin  by  the  exciting  wave. 


If  the  far  field  patterns  radiated  by  unit  dipole  moments  are 
denoted  by 
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where  the  superscripts  denote  the  type  and  direction  of  the  respective 
dipoles,  then  the  elements  of  the  polarizability  matrix  are  given  by 
the  reactions, 
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In  these  equations,  Eex,  etc.,  are  known  functions  of  (e,c)  on  the 
sphere  £  at  infinity  (from  reaction  matching)  sad  Eex,  etc.,  are  known 
functions  of  (0,<?>)  as  well  (from  Eq.  (6)),  so  it  is  only  a  matter  of 
numerically  integrating  their  dot  product  over  £  to  find  the  elements 
of  the  polarizability  matrix.  In  general  21  such  reaction  inteqrals 
must  be  evaluated. 


Typical  non-zero  values  of  the  polarizability  matrix  elements 
are,  for  the  examples  cited  in  the  text  of  this  report, 


Fig.  1: 
Fig.  2: 
Fig.  3: 


aez  =  1 .771xl0-4-j8.796xl0-7 
=  1.157xl0"3-j3.735xl0“5 


c iH  =-l-615xlO'2-j2.788xlO'2 


Fig.  4: 
Fig.  5: 


Fig. 

Fig. 


6: 

7: 


ex  ez  ex  ez 

x  ~  a  =  a  =  a  - 

ex  ex  ez  ez 

8 . 809x1 0_5-j4 . 374xl0-7 

ex  ez  ex  ez 

x  -a  =  a  -  a  - 

ex  ex  ez  ez 

5 ,666xl0”4-jl . 829x1 0”5 


Fig.  8:  a!*  =  a®y  =  7 . 056x1 0”4- j 1 .434x10 

Fig.  9:  fex  ey 

aZ  =  -3.077x1 0"4-j2 .666xl0-6 
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